Evolution of spectral function in a doped Mott insulator : surface vs. bulk 

contributions 
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We study the evolution of the spectral function with progressive hole doping in a Mott insulator, 
Lai-^Caa; VO3 with x = 0.0 - 0.5. The spectral features indicate a bulk-to-surface metal-insulator 
transition in this system. Doping dependent changes in the bulk electronic structure are shown to 
be incompatible with existing theoretical predictions. An empirical description based on the single 
parameter, U /W, is shown to describe consistently the spectral evolution. 
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It is well understood that the insulating ground state 
of transition metal compounds with integral occupancies 
of 3d levels (Mott insulators) arise from strong electron 
correlation effects. Recent years have seen a phenomenal 
resurgence of interest in studying the properties of such 
systems doped with charge carriers representing frac- 
tional 3d occupancies, leading to the discovery of many 
exotic properties such as the high temperature supercon- 
ductivity and colossal magnetoresistance Q|. A con- 
siderable amount of effort has been directed to modelling 
these systems theoretically, which is crucial for under- 
standing such properties. In recent times, there are very 
specific predictions about the evolution of the spectral 
functions in doped Mott insulators based on the Hubbard 
model in the limit of infinite dimension p]4 l. Consider- 
ably different ansatz, such as slave boson Jo . dynamical 
mean field theories || and exact calculations |5j for fi- 
nite systems, yield qualitatively similar results. These 
predictions are most useful, since the spectral function 
can be directly measured experimentally by photoemis- 
sion spectroscopy and thus, provide a very convenient 
testing ground for the suitability of the model. In or- 
der to investigate the spectral evolution as a function of 
doping we have studied Lai-^Ca^VOa for x — 0.0 - 0.5. 
LaVC>3 is a Mott insulator with the d 2 electronic con- 
figuration. Ca substitution dopes holes into the system 
continuously varying the electron count from 2 to 1.5 per 
V ion for the compositions studied. Transport and mag- 
netic measurements Q have shown that the system is 
an antiferromagnetic insulator for x < 0.2, but is Pauli 
paramagnetic metal for x > 0.2. Moreover, there is no 
significant structural modification suggesting nearly con- 
stant bandwidth (W), with a little (< 10%) change in 
Hubbard U across the compositions studied ||. Hence, 
the change in the electronic structure in this system is 
primarily driven by the changes in the doping level and 
thus, provide a suitable case to investigate the spectral 
evolution in a doped Mott insulator. We observe the sur- 
face electronic structure in this system to be considerably 
different from the bulk one and provide an algorithm to 
extract the true bulk electronic structure from the total 
spectra using photon energy dependence of the spectral 



functions. The extracted bulk spectra show that existing 
theories based on the Hubbard model are insufficient to 
describe the experimentally observed spectral functions 
for various dopings. We provide an empirical approach to 
describe the observed evolution of the spectral function 
and discuss the implications of our findings. 

Polycrystalline samples of Lai_ x Caj;V03 were prepared 
as reported elsewhere Q , from congruently molten states 
giving rise to large grains with strong intergrain bonding 
and about 1% of impurity phases in the grain bound- 
aries. The x-ray photoemission (XP) and ultraviolet pho- 
toemission (UP) spectroscopic measurements were car- 
ried out at liquid nitrogen temperature in a vacuum of 
4xl0~ 10 mbar with resolutions of 0.8 eV and 0.08 eV 
respectively. The cleanliness of the sample surface was 
maintained by periodical scraping with an alumina file 
in situ and was monitored by O Is and C Is spectral re- 
gions in XP and the characteristic region (9 - 12 eV bind- 
ing energy) in UP measurements. Reproducibility of the 
spectra with repeated scrapings was confirmed for each 
composition. We believe that the spectra are representa- 
tive of angle integrated spectra, since they are obtained 
over a wide (± 10°) acceptance angle and were always 
independent of sample rotation about the surface nor- 
mal. The bulk spectra were primarily extracted from the 
XP spectra, which by their nature average over a large 
part of the Brillouin zone. We implicitly assume that the 
incident photon energy-dependent cross-sections and the 
combined resolution do not seriously change the spec- 
tra over our narrow (3 eV) binding energy window; we 
substantiated this assumption explicitly for LaVC>3. The 
derived spectra for the bulk, within our resolution and 
over the relatively large energy range investigated here, 
are not affected by the disorder due to scraping. The bulk 
and surface compositions were found to be in agreement 
with nominal compositions from energy dispersive x-ray 
analysis and the intensity ratio of the core level spectra 
of the componenets, respectively. 

In Fig.l we show the He 1 UP spectra for all the com- 
pounds near the V 3d emission region. Each spectrum 
is dominated by a feature centered at about 1.5 eV be- 
low Ep , with no intensity at the Fermi level (Ep) for x 
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= 0.0 and 0.1. This feature is normally termed as in- 
coherent feature |l(J being the spectral signature of the 
lower Hubbard band (LHB) and corresponds to electron 
states essentially localized due to electron correlations. 
The spectral feature growing at Ep (coherent peak) with 
increasing x (x > 0.2) represents the delocalized conduc- 
tion electrons. Thus, these observations show a metal- 
insulator (MI) transition at x = 0.2 in agreement with 
the transport properties (8). However, the complete dom- 
inance of the incoherent feature suggests an overwhelm- 
ing presence of correlation effects nearly localizing charge 
carriers even in the x = 0.5 sample in contrast to the ob- 
served physical properties ||. As we shall show now, 
at least a part of this discrepancy arises from a surface 
transition of the electronic structure in these systems. 

The surface sensitivity of these electronic spectra can 
be varied by changing the exciting photon energy, since 
the escape depth (A) of the electrons depends sensitively 
on their kinetic energy We use this fact to delin- 

eate the surface and bulk electronic structures by record- 
ing the spectra using He i (21.2 eV), He n (40.8 eV) 
and monochromatized Al Ka (1486.6 eV) sources; among 
these the spectrum excited with He n radiation is ex- 
pected to be the most surface sensitive, while the Al Ka 
spectrum the most bulk sensitive. We show the spec- 
tra for each of the metallic compositions in Fig. 2, after 
subtracting the O p contributions (solid line in the inset 
of Fig.l) appearing at binding energies higher than 3.5 
eV following the procedures in Refs. jl2],[l3|]. For each 
sample, the spectrum recorded with He n radiation ex- 
hibits the single incoherent spectral signature with no 
intensity at Ep indicating absence of any coherent state. 
This clearly demonstrates that the electronic structure 
in the surface region of these samples remains localized. 
This most probably arises from an enhanced correlation 
effect near the surface region due to the reduced dimen- 
sionality at the surface and/or subtle changes of surface 
geometry (surface reconstruction) and consequently a re- 
duced bandwidth, even for x = 0.5 sample. This is not 
entirely surprising since the critical concentration even 
for the bulk MI transition in closely related, but slightly 
distorted, series Yi_ K CazV0 3 is indeed x = 0.5 @. The 
finite spectral weight in the coherent feature in the He i 
spectrum increases further in the Al Ka excited signal for 
each of the samples, consistent with an increasing bulk 
sensitivity of the technique. This substantiates further 
the view that the bulk electronic structure represents a 
metallic state at these compositions, though the surface 
layer has only localized states. Hence, to probe the true 
bulk spectral function, it is however necessary to separate 
the surface contribution from the total spectra in Fig. 2. 

The total normalized spectral intensity F(e) from the 
sample at normal emission is given by 

F(e) = F s (e)(l - e - d / x ) + F b {e) e - d/x 

where, F s (e) (given by the He n spectra) and F h (e) are 



the normalized surface and bulk spectral functions re- 
spectively and d is the thickness of the insulating surface 
layer. Expressing e d l x as a, we readily find that F b (e) is 
given by [a.F(e) — (a — l).F s (e)]. If a, which is a function 
of the photon energy via A, would be known, it would be 
straight forward to obtain the bulk spectrum, F b (e), from 
the above equation. In absence of a prior knowledge of 
a however, we note that the spectral shape of F b (e) rep- 
resenting the bulk electronic structure is independent of 
the photon energy. Thus, the extracted F b (e) from the Al 
Ka spectra must be consistent with that extracted from 
the He i spectrum. This consistency requirement can be 
cast into a least-squared-error procedure to estimate the 
values of a corresponding to He i and Al Ka photon en- 
ergies. Thus obtained values of and = 2.95 and axps 
= 1.38 correspond to d/Xxps = 0.32 and d/Xnei = 1-08 
leading to a reasonable Xxps/^Hei ~ 3.4. Moreover, 
using the estimates of A from the published literature 
[jll],[l5 16 1, the present result suggests that the thickness, 
d, of the insulating surface layer is approximately 8 A, 
which is roughly the dimension of the unit cell. 

In inset I of Fig. 3 we demonstrate the subtraction pro- 
cedure for obtaining the bulk spectra E^j. Thus ob- 
tained bulk electronic structures for the metallic compo- 
sitions are shown in Fig. 3 by open circles. The spectral 
weights show a systematic increase at Ep with increas- 
ing Ca, consistent with increasing conductivity in the 
series. However, the V 3d emission is seen to be spread 
to about 3 eV below Ep with the highest intensity in the 
vicinity of 1.5 eV. Our band structure calculations for 
both LaVC>3 and CaV03 however suggest that the elec- 
tronic density of states in the occupied part has a peak 
at Ep with a total spread of about 1 eV. Thus, it is clear 
that the bulk electronic structures as characterized by 
the spectra in Fig. 3 are strongly influenced by the elec- 
tron correlation effect, leading to substantial incoherent 
spectral signatures giving rise to the peak near 1.5 eV 
and an extension of the emission down to 3 eV below 
Ep. The doping-dependent spectral function has been 
calculated within the single band Hubbard model both 
by dynamical mean field approximations |(| and exact 
diagonalization techniques . These results indeed sug- 
gest the presence of the incoherent and coherent features 
in the spectra. However, any reasonable value of U/W 
with finite doping shows a much weaker contribution of 
the incoherent feature as compared to our experimental 
results. This weakened effect of correlation in the pres- 
ence of large doping within a single band Hubbard model 
is understandable, since such doping reduces the average 
number of electrons per site considerably from 1.0. In 
the present case, there are 1.5 electrons per V site even 
for the highest doped case (Lao.5Cao.5VO3). In order to 
understand the effects of correlation in such a multiband 
case, we have calculated the spectral function for a de- 
generate Hubbard model p8[ , since such results are not 
available in the literature. We show the resolution broad- 
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ened spectral functions for various values of U/W in inset 
II of Fig. 3 exhibiting a complete dominance of the co- 
herent feature and the presence of only weak incoherent 
features above 1.0 eV. Thus, the inability of the Hubbard 
model to account for the experimentally observed intense 
incoherent spectral features in these doped Mott insula- 
tors cannot be attributed to the multiband nature 
of the problem. 

In this context, it is important to realize that the end 
members of Lai-zCa^VOs have fundamentally different 
electronic properties, one (LaVOs) being a Mott insu- 
lator, while the other (CaVOs) is a correlated metal. 
Moreover, we note that the peak position near 1.5 eV 
and its extension down to 3 eV is remarkably similar to 
the spectral signature of LaVC>3 (see Fig. 3). The spec- 
tral weight shifted nearer to Ep in the metallic samples 
compared to LaVC>3 is clearly the signature of the co- 
herent feature growing in relative intensity with progres- 
sive doping of the Mott insulator. It turns out that we 
can model this coherent part of the spectral signature 
with ab initio band structure calculations for CaVC>3 
corrected for correlation effects within the perturbative 
treatment of the self-energy [p~S| ] . Considering linear con- 
tributions of the spectral features arising from LaV03 
and CaV03, the general sum rule for the conservation 
of electron number completely determines the relative 
contributions. This can be expressed as the total mea- 
sured spectrum p(e) = 2(1 — x)pLa(t) + xpca(f)- Here, 
Pia(e) is the normalized spectral function of the LaVC>3 
part and is completely determined from the experimen- 
tal XP spectrum of LaVC>3. x is also fixed by the com- 
position of the compound Lai_ x Ca x V03 and ensures a 
total electron count per V ion equal to (2-x). pca(e) 
is calculated from the band structure results corrected 
up to second order in U /W . Thus, p(e) depends on a 
single parameter, U/W, with all other contributions be- 
ing a-priori fixed. Noting the possible limitations of any 
parametric fitting procedures, we fix even the value of 
this parameter (U /W = 0.5) by requiring that the band 
structure results for CaV03 when corrected for correla- 
tion effects via self energy calculation, yields the correct 
m*/nib = 2.1 ^(j). Thus, the calculated total spectrum 
p(e) is completely devoid of any adjustable parameter. 
The change in the spectral features is determined solely 
by the value of x. The resulting p(e) are shown in Fig. 3 
by the solid lines for each of the compositions establish- 
ing a remarkable agreement between the experiment and 
the parameter-free simulation. This result suggests that 
LaVC>3 and CaVC>3 somehow retain their characteristic 
electronic structures even for the intermediate composi- 
tions indicating a dominance of local interactions in de- 
termining the electronic structures; this is however a rea- 
sonable expectation for strongly correlated systems. 

In conclusion, it is evident from our results that the sur- 
face electronic structure of early transition metal oxides 
can be qualitatively different from the bulk one. This 



realization is essential in order to critically discuss and 
evaluate the experimental electronic structure in terms 
of the existing many-body theories and various bulk sen- 
sitive low-energy properties. We present an algorithm to 
extract the bulk related spectra from the total spectra us- 
ing photon energy dependent measurements. It is shown 
that all known theoretical results based on the single 
band Hubbard model [||J7[] as well as the present multi- 
band Hubbard model, are clearly insufficient to explain 
the experimentally observed bulk spectral functions for 
any value of U /W. It is then remarkable and significant 
that a simple additive description of the spectral func- 
tions of the end members, LaVC>3 and CaVC>3, though 
empirical and speculative, is qualitatively and quantita- 
tively successful for all the compositions without using 
any adjustable parameter and therefore deserves serious 
considerations in terms of rigorous microscopic theories 
for such systems. A direct consequence of this obser- 
vation would be that the electronic structure of such a 
system is intrinsically inhomogeneous and cannot be de- 
scribed within a homogeneous model, such as the Hub- 
bard model. This is in contrast to the usual practice 
in the recent times to interpret the properties of doped 
Mott insulators in terms of the Hubbard model. While 
the presented empirical approach appears to be consis- 
tent with the spectral evolution as a function of doping 
of the Mott insulator including both the coherent and 
the incoherent parts, it remains to be seen if in future a 
proper theoretical model would provide a rigorous basis 
for understanding the physical properties of these sys- 
tems. 
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I. FIGURE CAPTIONS 

Fig.l He i spectra of Lai_ x Ca x V03. Expanded (x3) 
view of near Ep region is also shown. Inset shows ex- 
pected tailing of O p states. 

Fig. 2 Normalized (equal area) He i, He n and XP va- 
lence band spectra of Lai-xCa^VOs for x = 0.3, 0.4 and 
0.5. 

Fig. 3 Bulk spectra (open circles) of Lai-zCa^VC^. The 
solid lines show the synthesized spectra. The correspond- 
ing spectra of LaV0 3 (solid circles) and CaV0 3 (dashed 
lines) are also shown. Inset I shows the surface and bulk 
components of the spectral function for He i and XP spec- 
tra for x = 0.5(see jl7|]). Inset II shows the results ob- 
tained from multiband Hubbard model for various values 
of U/W. 
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